Dynamics of polymer translocation into a circular nanocontainer through a nanopore 
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Using Langevin dynamics simulations, we investigate the dynamics of polymer translocation into 
a circular nanocontainer through a nanopore under a driving force F. We observe that the translo- 
cation probability initially increases and then saturates with increasing F, independent of (f>, which 
is the average density of the whole chain in the nanocontainer. The translocation time distribution 
undergoes a transition from a Gaussian distribution to an asymmetric distribution with increasing 4>- 
Moreover, we find a nonuniversal scaling exponent of the translocation time as chain length, depend- 
ing on <j) and F. These results are interpreted by the conformation of the translocated chain in the 
nanocontainer and the time of an individual segment passing through the pore during translocation. 

PACS numbers; 87.15.A-, 87.15.H- 



I. INTRODUCTION 

Biopolymer translocation across a membrane through 
a nanopore has attracted broad interest because it is es- 
sential to many biological processes, such as DNA and 
RNA translocation across nuclear pores, protein trans- 
port through membrane channels, virus infection 
In a seminal paper, it has been experimentally demon- 
strated that an external applied electric field can effec- 
tively drive a single-stranded DNA and RNA molecules 
through the a-hemolysin channel and that the passage 
of each molecule is signaled by a blockade in the channel 
current Ql- Due to its numerous revolutionary applica- 
tions, including rapid DNA sequencing, filtration, gene 
therapy, fiovif-injection problems, and controlled drug de- 
livery [sl, m , considerable experimental [1, [sl-flTj , theoret- 
ical and numerical [T8l - l4l| studies have been devoted to 
polymer translocation. 

Based on both the basic physics as weW as a technology 
design perspective, the scaling of the average transloca- 
tion time T with the polymer length iV, t ~ iV", is an 
important measure in the study of polymer transloca- 
tion. The scaling exponent a reflects the translocation 
efficiency. Standard equilibrium Kramers analysis |42| of 
diffusion through an entropic barrier discovers r ^ N'^ 
for unbiased translocation and t ~ iV for driven translo- 
cation (assuming friction is independent of N) [ll, [l^ . 
However, the scaling behavior t ^ iV^ is incorrect at 
least for a self-avoiding polymer. The reason is that the 
translocation time is shorter than the Rouse equilibra- 
tion time of a self-avoiding polymer, th ~ N^^'^'^ , where 
the Flory exponent ly — 0.588 in three dimensions (3D) 
and 0.75 in two dimensions (2D) [i^ |4^. In contrast, 
the scaling behavior tr ^ N^^'^" for both phantom and 
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self-avoiding polymers was demonstrated by Chuang et 
al. [l^l and Luo et al. based on a 2D lattice model, 
the Fluctuating Bond and Langevin Dynamics (LD) [22j 
models with the bead-spring approach, respectively. For 
driven translocation by a transmembrane electric field, 
a lower bound r ^ N^^'^ for the translocation time 
was provided by Kantor and Kardar [2l|. Recent re- 
sults [2^, [2^ further show that there is a crossover from 
T ~ N^-^"^ for faster translocation processes to t ^ N^^'^ 
for slower translocation in 3D. 

In addition to a transmembrane electric field, the driv- 
ing force can also be provided by a pulling force exerted 
on the end of a polymer [21], [4^ , binding particles (chap- 
erones) j46l . l47j , longitudinal drag flow in a nanochannel 
[isj , or geometrical confinement of the polymer (49l - [53| . 

For driven translocation by a transmembrane elec- 
tric field, however, above physical pictures are based on 
translocation into an unconfined trans side. Very little 
attention is paid to the dynamics of translocation into 
confined environments (ssj . Even for the case of an un- 
confined trans side, the translocated chain is highly com- 
pressed during the translocation process because it does 
not have enough time to diffusion away from the exit 
of the pore under fast translocation conditions [H, [2^ , 
and thus even more severe nonequilibrium effects are ex- 
pected under confinement. 

Polymer translocation into a unclosed confined envi- 
ronment, such as a slit or a fluidic channel [s^, is sub- 
ject to a large entropic penalty due to reduced acces- 
sible degrees of freedom, which should dramatically af- 
fect the translocation dynamics. In addition, studies 
on translocation into conflned geometries will shed light 
on the dynamics of the packaging of DNA inside virus 
capsid |5J-[5g|, where the entropic penalty is from both 
the limited volume and the additional self-exclusion of 
the chain. Using the stochastic rotation dynamics (also 
called multiparticle collision dynamics) simulations [STj . 
Ali et al. [Sg have investigated the packaging of flexible 
and semiflexible polymers. However, they only focused 



2 



10 


• wall particle 
polymsr particlB 


(a) 


> 


€fh f 




-10 







-40 -30 -20 ^ -10 



>-0 



• wall particle 
o polymer particle 


(b) 


10 


• wall particle 
o polymer particle 


(c) 






>- 










-10 







-40 -30 -20 -ID ID -40 -30 -20 -10 

X X 



FIG. 1: Schematic representation of polymer translocation 
through a pore into a 2D spherical compartment with radius 
R through a nanopore driven by an external electric force 
F = 5: (a) before translocation, (b) during translocation, and 
(c) after translocation. The width of the pore is w = 1.6a. 



on the packing rate and the pauses in the packaging pro- 
cess. Many important details are still not clear, such as 
(a) the translocation probability, (b) the distribution of 
the translocation time, and (c) the average translocation 
time as a function of the density of the chain, the chain 
length, and the driving force. To this end, we investigate 
the dynamics of polymer translocation into 2D circular 
nanocontainer using Langevin dynamics simulations. 

The paper is organized as follows. In Sec. II, we briefly 
describe our model and the simulation technique. In Sec. 
Ill, we present our results. Finally, the conclusions and 
discussion arc presented in Sec. IV. 



II. MODEL AND METHODS 

In our simulation, we model the polymer chains as 
bead-spring chains of Lennard- Jones (LJ) particles with 
the Finite Extension Nonlinear Elastic (FENE) poten- 
tial [13. Excluded volume interactions between beads is 
modeled by a short range repulsive LJ potential: 



nanocontainer through a nanopore driven by a trans- 
membrane electric field. The circular nanocontainer with 
a pore of width w ~ 1.6cr is formed by stationary wall 
particles within a distance a from each other. Between 
all monomer- wall particle pairs, there exits the same 
short range repulsive LJ interaction as depicted above. 
In Langevin dynamics simulation, each monomer is sub- 
jected to conservative, frictional and random forces, re- 
spectively. 
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where m is the monomer's mass, ^ is the friction co- 
efhcicnt, is the bead's velocity, and Ff is the ran- 
dom force satisfying the fluctuation-dissipation theorem. 
The external force is expressed as 'Fext = Fx, where F 
is the external force strength exerted exclusively on the 
monomer in the pore, and i is a unit vector in the direc- 
tion along the pore axis. This driving force mimics the 
role of the portal molecular motor for packaging of DNA 
molecules into capsid. 

In the present work, we use the LJ parameters e, ct, 
and 771 to fit the system energy, length and mass scales, 
respectively. As a result, the corresponding time scale 
and force scale are given by t^j ~ (mcr^/e)^/^ and e/a 
respectively, which are order of ps and pN, respectively. 
We set ksT ~ 1.2e, which means the interaction strength 
e to be 3.39 x 10^^^ J at actual temperature 295 K. This 
leads to a force scale of 3.3 pN for a chain with segmental 
length (7 = 1 nm. In our simulation, the dimensionless 
parameters are chosen to be i?o = 2, A: = 7, ^ = 0.7. 
Then, the Langevin equation is integrated in time by the 
method proposed by Ermak and Buckholz [60j . 

Initially, the first monomer is placed at the pore center, 
as shown in Fig. [IJa), while the remaining monomers are 
undergoing thermal collisions described by the Langevin 
thermostat to obtain an equilibrium configuration. The 
translocation time t is defined as the time duration be- 
tween the beginning of the translocation and the last 
monomer of the chain entering into the circular nanocon- 
tainer, see Fig. [Ijc). Typically, the averaging is done over 
2000 successful translocation events. 
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for r < 2I/V and for r > Here, a is the di- 

ameter of a monomer and e is the potential depth. The 
connectivity between neighboring beads is modeled as a 
FENE spring with 
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where r is the distance between consecutive monomers, 
k is the spring constant and Rq is the maximum allowed 
separation between connected monomers. 

As shown in Fig. [U we consider a two-dimensional 
(2D) geometry, where the chain is translocating into a 



III. RESULTS AND DISCUSSIONS 
A. Theory 

The conformation of a linear chain confined in a space 
of characteristic size lower than its radius of gyration 
has been broadly studied [6ll - [66j . The scaling behav- 
ior of free energy for polymer in full three-dimensional 
confinement, such as a spherical cavity, is different com- 
pared with slit-like confinement and tube-like confine- 
ment [6l,[63,[6^. 

Following Refs. [6^and[66l we first consider the equilib- 
rium behavior of a linear chain of length confined into 
a two-dimensional circular nanocontainer of the radius 
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R ^ Rg. Here, Rg '-^ N'^'^^a is the radius of gyration 
of the chain without confinement, with = 0-75 be- 
ing the Flory exponent in 2D [i^, ll^l- The notation 
means that two quantities are proportional to each other. 
The density of the bead in the nanocontainer is 

<j>^ Na^/{2Rf. (4) 

According to the blob picture, the chain will form blobs 
with size ^6, and each blob contains g ~ (^b/f)^/''^" 
monomers. The density of the monomer in a blob is 
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which should be the same as that of the bead in the whole 
nanocontainer, (f>biob = 0- Therefore, based on Eqs. (jj]) 
and (O, the blob size is 

6 cr0-''^o/(2''2"-l) = C^(/.-l■^ (6) 

which indicates that the blob size decreases with increas- 
ing (j>. 



Then, the number of blobs is 



Ub = N/g - Ar0i/(2^2D-i) ^ 



(7) 



The free energy cost is proportional to the number of 
blobs, and thus in units of /cgT it is 
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This shows that the free energy of the chain increases 
with 0^. 

For polymer translocation into a circular nanocon- 
tainer, an external driving force F is necessary to over- 
come the entropic force from the already translocated 
beads at cis side. With increasing the time, the den- 
sity of the chain in the nanocontainer increases. Once 
the translocated portion of the chain feels to be con- 
fined, it will form blobs and the blob size becomes smaller 
and smaller with the time. Thus, the entropic repulsion 
f{4>{t)) exerted by the already translocated monomers 
increases with the time during the translocation process. 
Owing to the highly nonequilibrium nature of the translo- 
cation process, it is difficult to calculate f{<f>{t)). 

To examine the the effect of the confinement on the 
translocation time, we define (f> as the chain density in the 
circular nanocontainer after the translocation and 
as the average resisting force for the whole translocation 
process. Then, we can write the translocation time as 

T ^ Ny[F - ^ Ny[F{l - (9) 

where F is the external driving force in the pore and a 
is the scaling exponent of r with chain length N. Here, 
F{1 ~ f{(f>)/F) is the effective driving force. For an un- 
confined translocation (R — >■ oo), translocation time is 

Too ~ iV"/F, (10) 

and then we obtain 
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Based on this relationship, we can explore the depen- 
dence of the average resisting force f{(p) on the density 
of the chain 6. 
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FIG. 2: (a) Translocation probability as function of the driv- 
ing force F for A'' = 32, 128 and different density of the chain 
(f). (b) Translocation probability as function of the density of 
the chain 4> for chain length iV = 128 {F = 2, 3, 5, 10 and 15) 
and Ji = 9.5 (F = 2). 



B. Numerical results 

During the translocation process, the effective driving 
force is F{1 — f {(/))/ F). In our simulations, if F < 1, the 
translocation probability is too small and the transloca- 
tion time is too long, particularity for higher 0. To obtain 
effective statistical average for translocation probability 
and translocation time, we present results for F > I. 



1. Translocation probability as a function of the driving 
force F and the density of the chain (f> 

We refer to a successful translocation as the event when 
the chain fully enters into the circular nanocontainer, as 
shown in Fig. [He). Otherwise, the translocation is con- 
sidered as unsuccessful. Fig. [H^a) depicts the transloca- 
tion probability, denoted as Ptrans , as a function of the 
driving force F for different chain lengths {N = 32 and 
128) and density of the chains (0 = 0.15, 0.3 and 0.45 
through the change of R). With increasing the driving 
force _F, Ptrans iucrcascs rapidly first, and then slowly 
approaches saturation at larger F, almost independent 
on N and (f). As shown in Fig. [IJb), we also observe 
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FIG. 3: Number of beads packaged against time for the chain 
length A'' = 128, the driving force F = 2, and different <j}. 
Here, all the data points are not averaged, and are from a 
typical successful translocation event. The pause is defined 
as the situation that the number of packaged beads is constant 
with time, and the typical pause is marked by arrow in the 
figures. The average time duration of pauses for cp = 0.45 is 
longer than those for (j> — 0.2. 



that for a fixed N value, Ptrans is independent of we 
examined for different driving forces. 

If the radius of the circular nanocontainer R is fixed, 
(j) changes with N . Intuitively, increasing the density of 
the chain will lead to the decrease of Ptrans due to the in- 
crease of the resisting force. However, we find that Ptrans 
is also independent of for R = 9.5 and F = 2. Dur- 
ing the translocation process, the density of the chain (jj 
of translocated monomers in the nanocontainer increases 
with the time. But the effect of (p at the early time of 
the translocation process is negligible, which dominates 

Ptrans • 
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FIG. 4: Histogram of translocation times for the chain length 
— 128 and different F and 4>. Here, F — 2 and 5 and 4> = 
0.2 and 0.45. The inset shows the distribution of translocation 
times for F = 2 and <6 = 0.45. 



2. Distribution of the translocation time 

Fig. |4] shows the distribution of translocation times for 
N — 128 and different F and 0. Under the driving force 
F = 5, the histograms for both (j> = 0.2 and 0.45 show 
narrow Gaussian distributions, but the position of the 
peak moves to a longer time for = 0.45 compared with 
that for (j> = 0.2. With decreasing the driving force to 
F = 2, for (/) = 0.2 the histogram still looks like Gaussian 
distribution but is wider by comparison with the cases 
for F = 5. However, for cj) = 0.45 the distribution is 
much wider and highly asymmetric, see the inset of Fig. 
m These distributions can be understood by taking into 
account the effective driving force F{1 ~ f{(p)/F). With 
increasing </>, the effective driving force F{1 ~ f{(j>)/F) 
decreases, which leads to the peaks of the distributions 
moving to longer times. For F = 2 and (j) = 0.45, the 
effective driving force becomes quite small during the late 
time of the translocation process, resulting in an highly 
asymmetric distribution. 



Even if (jj is large, which means the resisting force 
is strong at the late time of the translocation process, 
the chain can experience backward motion to decrease 
the resisting force. Thus, it is still difficult for the first 
monomer to exit the pore, and Ptrans is independent of 
(j>. This is can be seen in Fig. [31 where a typical plot of 
the number of beads packaged against time for a success- 
ful translocation event. For the chain length N = 128 
and the driving force F = 2, we observe pauses for both 
4> = 0.2 and 0.45, where the number of packaged beads is 
constant with time. The average time duration of pauses 
for (/) = 0.45 is longer than those for (j> = 0.2. We should 
point out that pauses are also observed in the stochas- 
tic rotation dynamics simulations |58l|. as well as in the 
experiments on DNA packaging |54| . 



3. Translocation time as a function of the density of the 
chain (j) 

Fig. EJa) shows the translocation time r as a function 
of the density of the chain (f> under the driving force F = 2 
for a fixed chain length = 128 or a fixed radius of 
the nanocontainer R = 9.5. For a fixed chain length 

= 128, we increase by decreasing R. For this case, 
we find that r initially increases very slowly with 0, and it 
increases very rapidly after a critical density of the chain, 
(j)c ~ 0.4. For a fixed R = 9.5, 4>\s increased by increasing 
N. For (j) < 0.2, we find r - ^^-^^ implying r - N^-'^'^. 
The scaling exponent is close to 21/ for translocation into 
unconfincd system, reflecting that cj) does not affect the 
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FIG. 5: (a) Translocation time r as a function of the density 
of the chain, 0, under driving force F = 2 for a fixed chain 
length 128 or a fixed radius of the nanocontainer R = 9.5. (b) 
1 — (too/t) as a function of (j>. 



scaling exponent for cj) < 0.2 due to F{1 — f {(/))/ F) w F. 
However, r rapidly increases with cj) particularly after (f>c 
because f{4>) becomes important. 



Fig. [5{b) shows 1 — {too/t) as a function of (p for F = 2. 
Here, (f> varies by either changing R or N. The data 
points for both cases almost collapse on one curve. With 
increasing (/), 1 — (too/t), which is proportional to f{(j))/F, 
increases slowly at first, followed by a rapid increases. For 
strong confinement = 0.45, the average resisting force 
/((/)) is about 80% of the driving force F. This indicates 
that the effective driving force F{1 — f{(j))/F) is quite 
small (2 X 20% = 0.4), leading to asymmetric distribution 
with long tail of translocation times shown in the inset of 
Fig. H From the inset of Fig. [S^b), - c/)^ , where 

P is very small for the weak confinement {(j) < 0.2), 2.5 
for the moderate confinement (0.2 < cf) < 0.4) and 7.1 
for the strong confinement {(j) > 0.4), reflecting the very 
rapid increase of the resisting force with (f> for larger 0. 
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FIG. 6: Translocation time r as a function of the chain length 
A*' under the driving force F = 2: (a) for different density of 
the chain, ij>, (by changing _R), and (b) for R = 9.5. 
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FIG. 7: Translocation time t as a function of the chain length 
for different density of the chain (j> and different chain 
lengths A'' under the driving force F = 5. 



4- Translocation time as a function of the chain length N 

Fig. [6l^a) shows r as a function of the chain length 
under the driving force = 2 for different density of 
the chain 0. In the simulations, for a fixed (f> we changed 
R for different A''. We find r ^ N" with a being the 
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scaling exponent. For a weak confinement, <p = 0.1 we 
observe a = 1.47, in agreement with the previous study 
for translocation into an unconfined system p2j . where 
a = 2zy = 1.5 for TV < 200 in 2D. This resuh also in- 
dicates that the weak confinement does not change a 
due to the effective driving force F{\ — f{(j>)/F) w F. 
However, with increasing (/), the dynamics is different for 
translocation into a nano container. We observe that a 
slightly decreases first from a ~ 2i/ and then increases 
with increasing </), such as a = 1.65 for 4> = 0.45. With 
increasing </) from 0.40 to 0.45, the resisting force /(0) 
increases rapidly with 0, leading to a great decrease of 
the effective driving force F{1 — f{(f>)/F). This is the 
reason for the increase of a. Here, the exponent a in- 
creases with decreasing the effective driving force, which 
also confirms previous results on forced translocation in- 
side an unconfined region under different driving forces 
j29l \q1\ and results on translocation into two-dimensional 
fluidic nanochannels for a fixed driving force and differ- 
ent channel widths [s^. Recently, Lehtola et al. [gl] 
also observe a nonuniversal scaling exponent for forced 
translocation inside an unconfined region, however, the 
exponent a increases with increasing the driving force. 

For a stronger driving force = 5, the scaling expo- 
nent a almost docs not change even for = 0.45, see 
Fig. [71 This is because the driving force is dominant 
compared with the resisting force during the transloca- 
tion process, namely, F{\ — f{(j>)/F) F is still correct 
for larger F. 

As a comparison. Fig. El^b) shows r as a function of 
the chain length N under the driving force F = 2 for 
a fixed radius of the nanocontainer R — 9.5, where the 
resisting force f{(j)) increases slowly with for short N, 
and thus we still observe a « 2i/ for short chains. For 
longer chains, the resisting force /(0) increases rapidly 
with TV due to both the limited volume and the additional 
self-exclusion of the chain, leading to a much more rapid 
increase of the translocation time t. 



5. Translocation time as a function of the driving force F 

In previous study for translocation into the unconfined 
system [13, [1^, it has been shown that t ~ F~^ for lower 
forces and r F~^-^ for higher driving forces. 

For translocation into confined nanocontainer, the 
chain density in the nanocontainer, 0, increases with 
time, leading to that the resisting force f{4>) increases 
with time for a typical translocation event. In addition, 
the value of /(0) increases with cj) as shown in Fig. [Hb). 
For the lower density of the chain such as cj) = 0.1, the 
effective driving force F{1 ~ /{(/)) /F) « F, and thus we 
observe the similar results, see Fig. [51 c ompared with 
translocation into unconfined systems [22|. For F < 5, 
we observe r ~ F~^ before crossing over to r ^ ^-0.76 
for F > 5. The scaling of t ^ F~^ can be under- 
stood by considering the balance of the frictional damp- 
ing force with the driving force F. For F > 5, a dif- 
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FIG. 8: Translocation time r as a function of the driving force 
F for chain lengths N = 128 under the density of the chain 
= 0.1 and 0.45 in 2D. 

ferent scaling behavior is observed due to a pronounced 
non-equilibrium effect. However, for = 0.45, r rapidly 
decreases with F for F < 4 because f{(f>) is important 
and cannot be negligible for the translocation dynamics. 
For F > 5, f{(f>) becomes less and less important due to 
F(l - f{(t))/F) K. F for larger F, and thus r - F^^ is 
observed. 

Fig. [9l shows the radius of gyration of the chain, Rg , 
at the moment just after the translocation for N = 128 
as a function of the density of the chain (j) under different 
driving forces. For the weak confinement regime with 
small density of the chain 0, the radius of gyration of 
the translocated chain Rg under the strong driving force 
is much smaller than that under the weak driving force. 
The reason is that the translocated chain are more com- 
pressed nearby the entrance for strong driving forces, see 
the chain conformations for N ~ 256 and = 0.1 under 
F = 2 and F = 15 shown in Fig. [TUl 

For equilibrated chains, the radius of gyration of the 
chain in sphere is proportional to the radius of sphere 
under the strong confinement, namely, Rg ~ (N/cj))^^'^ ^ 
A^/^0~^/^. The numerical exponent of Rg with in- 
creases from -0.42 for F = 2 to -0.37 for F = 15. Due to 
the non-equilibrium effect, the chain is not homogenous 
distributed in the nanocontainer, and thus the simulation 
numerical exponent does not agree with -1/2. Particu- 
larly, for larger F, the segments of the chain are locally 
oriented along the direction perpendicular to the axis of 
the pore. 

6. Waiting time distributions 

It is necessary to explore the dynamics of a single seg- 
ment passing through the pore into the nanocontainer. 
The time of an individual segment passing through the 
pore is significantly affected by the non-equilibrium na- 
ture of the translocation process. We numerically calcu- 
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FIG. 9: The radius of gyration of the chain Rg at the mo- 
ment just after the translocation for N = 128 under different 
driving F force versus the density of the chain (p. 
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FIG. 11: Waiting time distribution for different densities of 
the chain (j). Here, the chain length is A*' = 128, and the 
driving force is _F = 2. 




FIG. 10: The chain conformation at the moment just after the 

translocation for the chain length iV = 256 and the density ^IG. 12: Waiting time distribution under different driving 
of the chain <^ = 0.1 under the driving forces f = 2 and 15 in ^^^'^''^ ^- The chain length N = 128 and density of the chain 



lated the waiting times for each monomer of the chain. 
The waiting time of bead s is defined as the average time 
between the events that bead s and bead s + 1 first pass 
through the exit of the pore. 

Fig. [TT] shows the waiting time distribution for the 
chain length N = 128, the driving force F = 2 and dif- 
ferent densities of the chain (j). The obvious feature is 
that the waiting times almost experience the same path- 
way for smaller s, due to lower densities of the chain in 
the nanocontaincr. As observed for translocation into the 
unconfincd system [22|, [2^ [3^ , we observe almost sym- 
metric distributions with respect to the middle monomer 
s = N/2 for weak confinements such as = 0.1, due to 
the effective driving force F{1 — f{4>)/F) w F. With in- 
creasing (/), the waiting times increase only after a certain 
s, which is a little different compared with previous stud- 
ies (ssj for polymer translocation into a two-dimensional 



fiuidic channel with a narrow channel width where the 
waiting times always increase. Particularly, it takes much 
longer time for beads s > N/2 to exit the pore. For 
(j) ~ 0.35, the waiting times approximately saturate after 
s > N/2. But for = 0.40 and 0.45, the waiting times 
increases very rapidly after s « 100 because /(</)) is im- 
portant and the effective driving forces F{1 — f{(j))/F) 
are greatly reduced. A big difference towards the end of 
the packaging is due to less backward motion for several 
end monomers, which leads to shorter waiting times. 

For strong driving forces, the effect of /(0) on the wait- 
ing times is reduced because even if = 0.45 the effective 
driving forces F{1 — f{(f>)/F) are more and more close to 
F with increasing F, see Fig. [T^ For F = 15, the waiting 
times arc almost not changed ioi cj) — 0.2 and 0.45. 
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IV. CONCLUSIONS 

Using Langevin dynamics simulations, we investigate 
the dynamics of polymer translocation into a circular 
nanocontainer through a nanopore under a driving force 
F. We observe that the translocation probability ini- 
tially increases and then saturates with increasing F, 
independent of (f), which is the average density of the 
whole chain in the nanocontainer. The translocation time 
distribution undergoes a transition from a Gaussian dis- 
tribution to an asymmetric distribution with increasing 
(j). Moreover, we find a nonuniversal scaling exponent of 
the translocation time as chain length, depending on (j) 
and F. These results arc interpreted by the translocated 
chain conformation in the nanocontainer and the time of 
an individual segment passing through the pore during 



translocation. 

Our findings have also shed light on dynamics of the 
packaging of DNA inside virus capsid. Certainly, it is still 
difficult to directly compare our results with DNA pack- 
aging because DNA is a semifiexible polymer which can 
get entangled and knotted in three dimensions. In future 
studies, it would be interesting to consider the infiuences 
of the shape of the nanocontainer and the stiffness of the 
chain on the translocation dynamics. 
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